Universality of the Kondo Effect in a Quantum Dot out of Equilibrium 
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We study the Kondo effect in a quantum dot driven out of equilibrium by an external ac field. 
The Kondo effect can be probed by measuring the dc current induced by an auxiliary dc bias Vdc 
applied across the dot. In the absence of ac perturbation, the corresponding differential conductance 
G(Vdc) is known to exhibit a sharp peak at Vdc = 0, which is the manifestation of the Kondo effect. 
In the equilibrium, there exists only one energy scale, the Kondo temperature Tk, which controls all 
the low-energy physics of the system; G is some universal function of eVdc/Tk- We demonstrate that 
the dot out of equilibrium is also characterized by a universal behavior: conductance G depends on 
the ac field only through two dimensionless parameters, which are the frequency u and the amplitude 
of the ac perturbation, both divided by Tk ■ We find analytically the large- and small- frequency 
asymptotes of the universal dependence of G on these parameters. The obtained results allow us to 
predict the behavior of the conductance in the crossover regime Twj ~ Tk- 

PACS numbers: 73.23.Hk, 85.30. Vw, 72.15.Qm 



I. INTRODUCTION 

The Kondo effect results from exchange interaction of 
itinerant electrons with a localized spin state. This inter- 
action leads to local spin polarization of the electron gas. 
The polarization becomes significant only at low, tem- 
peratures, due to the existence of collective statestl with 
small excitation energies. Simultaneously with the mod- 
ification of the spin susceptibility, the scattering prop- 
erties are modified significantly from those of the orig- 
inal localized spin state. The modification is especially 
striking in the case of antifcrromagnetic exchange inter- 
action, when the "spin cloud" formed out of free electron 
gas tends to screen the magnetic moment of the localized 
state. In this case, the scattering cross-section grows as 
the temperature is being lowered, and reaches the uni- 
tary limit at T — > 0. This phenomenon is responsible for 
the non-monotonic temperature behavior of the resistiv- 
ity of metals with magnetic impurities at low tempera- 
tures, which was the first experimentally observed mani- 
festation of the Kondo resonance.El However, the system 
of magnetic impurities embedded in a metal sample does 
not offer much control over the parameters even at the 
stage of the fabrication of the sample, not speaking of 
the dynamic variation of the parameters in the course of 
an experiment. 

Another class of systems, whose transport properties 
can also be affected by the Kondo effect, but which offer 
a much higher level of control over the system parame- 
ters, is provided by quantum dots. A quantum dot in a 
semiconductor planar heterostructure is a confined few- 
electron system contacted by sheets of two-dimensional 
electron gas (leads). If the total number of electrons on 
the dot is odd, then the dot is similar to a magnetic im- 
purity. Junctions between the dot and the leads produce 
an overlap of the states in the dot and in the 2D leads. 
This overlap leads to the exchange interaction between 
the spin of the dot and spins of the itinerant 2D elec- 



trons. At sufficiently low temperatures, the "spin cloud" 
is formed of the electrons in the leads. Like in bulk met- 
als, the scattering off the resulting many-body state is 
being enhanced as the temperature is being lowered, and 
reaches unitary limit at T — > 0. Different are only the 
manifestations of the Kondo resonance in bulk metals 
and in the quantum dot systems. In bulk metals, the 
enhancement of scattering by the Kondo resonance in- 
creases the resistivity. In quantum dot systems, on the 
contrary, the scattering facilitates transport through the 
dot. Kondo effect results in a specific temperature depen- 
dence of the linear conductance across the dot. If tem- 
perature is lowered, the conductance first drops due to 
the Coulomb blockade phenomeaonp and then increases 
again due to the Kondo effect.El At T — ► 0, the conduc- 
tance G reaches its maximum, which corresponds to the 
unitary limit of tunneling. Recently, the Kondo effect in 
a quantum dot was observed experimentally.!] □ 

The quantum dot devices are highly controllable, and 
can be operated in regimes inaccessible in the conven- 
tional magnetic impurity systems, that were used pre- 
viously for studying the Kondo effect. Irradiation of a 
quantum dot with an ac field offers a new, clever way 
of affecting its dynamics, which enables one to study 
the Kondo anomaly in essentially non-equilibrium con- 
ditions. The ac field can be applied to the gate, thus 
modulating the dot's potential with respect to the leads; 
alternatively, one may apply ac bias to the leads. In any 
case, driving the system out of equilibrium affects the dc 
conductance discussed above. Measuring the dc I — V 
characteristics, one can investigate the effects of the ir- 
radiation on the Kondo anomaly. 

A generic theoretical description of a quantum dot op- 
erates with a significant number of parameters and en- 
ergy scales describing the system. Nevertheless, in the 
case of no ac field, it turns out that the low-energy prop- 
erties of the quantum dot system which are related to 
the Kondo effect are controlled by only one relevant en- 
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ergy scale, which is the Kondo temperature Tk- The 
Kondo temperature, in turn, depends on the microscopic 
parameters of the system, e.g., on the gate voltage and 
conductances of the dot-lead junctions. Such universal- 
ity allows for easier understanding and description of the 
problem. 

The ac field introduces new parameters to the prob- 
lem, thus apparently breaking the universal description, 
which is valid in the static case. This re-emerging abun- 
dance of parameters makes it difficult to build a con- 
sistent description of the effects of the ac field on the 
Kondo conductance. The theoretical work performed up 
to the date, was concentrated on some specific regimes. 
Goldin and AvishaiQ considered the case of very strong 
ac bias with the help of the third-order perturbation the- 
ory in the dot-lead coupling. Nordlander et ala analyzed 
the effects of ac field of sufficiently high frequency, which 
ionizes the dot. They have conjectured that even at the 
temperature of the thermal bath T = 0, the finite rate 
of the dot ionization results in a finite effective temper- 
ature "seen" by the Kondo state. This way, irradiation 
provides the cut-off for the Kondo singularity and n 
duces the conductance. Later, we have demonstrated! 
that even in the absence of the dot ionization, irradi- 
ation is able to flip the spin of the dot, thus bringing 
decoherence into the Kondo state and diminishing the 
Kondo effect. In addition to analytical methods, a num- 
ber of numerical approaches have been usecol3 to study 
the conductance of a Kondo system out of equilibrium 
at certain sets of values of the bare parameters of the 
system. Because of the large number of the parameters 
involved, the results of such calculations are hard to ana- 
lyze. At any rate, such a limited consideration could not 
reveal the universalities of the problem. In our view, it 
also does not provide an insight into the regimes which 
do not allow for a perturbative treatment. 

In our earlier papeiiJ we have learned how to apply 
the Renormalization Group (RG) technique for a Kondo 
system out of equilibrium. It allowed us to sum up the 
infinite series of the perturbation theory in the dot lead 
coupling. This treatment, valid at relatively high deco- 
herence rates, yielded the expression for the conductance 
as a function of only one parameter, which is the ra- 
tio of the decoherence rate %/t to the Kondo tempera- 
ture. Thus we established that the Kondo temperature 
Tk remains a meaningful parameter of the theory, even if 
the ac field strongly suppresses the manifestations of the 
Kondo effect. We do not see a way, however, of measuring 
the decoherence rate directly. Moreover, the definition of 
this quantity in Ref. || makes sense only at a sufficiently 
high frequency of the ac field, Tuo ^> Tk- Therefore, the 
dimensionless ratio Ti/tTk can not be the only parameter 
describing the effect of ac field on the Kondo system. 

In this paper, we find the correct dimensionless vari- 
ables that characterize the amplitude and frequency of 
the ac field perturbing the Kondo system. If the ac field 
biases the dot, then the proper variables are eVac/Tif 
and Tiuj/Tk, where V ac is the amplitude of the ac bias. 



In the case of ac modulation of the gate voltage, the per- 
turbation introduced by the ac field is characterized by 
the corresponding dimensionless variation of the Kondo 
temperature 5Tk /Tk, which in principle can be indepen- 
dently measured. In terms of the proper pair of variables, 
the behavior of the Kondo conductance is universal. We 
find analytically the asymptotes of the universal depen- 
dence by further developing the RG treatment (valid in 
the case of strong suppression of the Kondo effect! and 
by generalizing the Nozieres' Fermi-liquid theoryliil onto 
the non-equilibrium case (which adequately describes the 
limit of weak perturbation by a low- frequency field) . 

Within this picture, we are able to describe in a con- 
sistent way the effect of irradiation in a wide range of 
frequencies of the ac perturbation - from zero to the dot 
ionization threshold; this includes the most interesting 
for the current experimentstj region hui ~ Tk- The 
remarkable feature of the presented description is that 
the Kondo temperature remains the only relevant energy 
scale, despite the essentially non-equilibrium character of 
the problem. 

The outline of our paper is as follows. In Sec. |l| we 
introduce the description of the system by means of the 
time-dependent Kondo Hamiltonian. Then we derive the 
expression for the Kondo conductance of the dot in the 
absence of ac field. The purpose of this derivation is to 
present the formalism which later would be suitable to 
describe the non-equilibrium states produces by the irra- 
diation. 



In Sections III- VI we consider the effect of ac mod- 



ulation of the gate voltage on the Kondo conductance. 
At higher frequencies, Tllj/Tk ^ 1, the decoherence of 
the Kondo many-body state is the principal channel via 
which the ac field influences the Kondo anomaly. This 
frequency domain is considered in Sections III-1V. De- 
pending on the strength of modulation, 5Tk/Tk, the 
suppression of the Kondo conductance is s ignificant (Sec- 
tion III), or relatively weak (Sections |lV|). In Sec. |v| we 
consider the limit of very small frequencies of the ac field. 
The decoherence probability in these conditions is expo- 
nentially small. However, the adiabatic evolution of the 
collective state, caused by the ac modulation, affects its 
scattering properties. It results in small deviation of the 
conductance from the unitary limit. The results of Sec- 
III V| are summarized in Sec. VI. 



tions 



In Sec. VII we consider the effect of ac bias on the 



Kondo conductance. It turns out that a strong suppres- 
sion of the Kondo effect is possible only if eV ac /TK ^> 1. 
Upon the increase of the frequency, the suppression di- 
minishes. This frequency dependence is opposite to the 
one in the case of gate voltage modulation. 



In Sec. |VIII| we consider the satellite peaks which can 
be created by irradiation in the I—V characteristic of a 
quantum dot. 

Finally in Sec. [X we compare experimental resultsEl 
with our theory. 
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II. THE KONDO EFFECT IN A QUANTUM DOT 

A. The model 

The system we study is a quantum dot attached to two 
leads by high-resistance junctions, so that the charge of 
the dot is nearly quantized. The Kondo effect emerges in 
a quantum dot occupied by an odd number of electrons 
at temperatures below the mean level spacing in the dot. 
Under such conditions, the topmost occupied level is spe- 
cial, since it is filled by only one electron. It is this level 
which produces the Kondo effect. The other levels, oc- 
cupied by two or no electrons are unimportant in our 
discussion (similarly to the inner shells of a magnetic im- 
purity in the conventional Kondo effect). Therefore the 
model of the dot attached to two leads can be truncated 
to the Anderson single-level impurity model, 

h = E & + ey «) C L* C fe<xa + E ( U « C LA + H - c -) 

k,cr,a k,<r,a 

+ y2(-E d + eV^ot cos wt)d£d ff + Ud\d^d\d^ ; (la) 

a 

V l .r = ±i [V dc + V ac cos(w't + <fo)] , 



r a = 2irv\v a \ 



(lb) 



Here the first two terms correspond to non-interacting 
electrons in the two leads (a = L,R), and tunneling of 
free electrons between the dot and the leads, respectively. 
The dot is described by the third and fourth terms of the 
Hamiltonian, Ed and U — Ed are the ionization and the 
electron addition energy, respectively. The tunneling ma- 
trix elements v a are related to the widths r a by Eq. (|lb|), 
where v is the density of states in a lead. The ac field 
can be applied to the gate, which is coupled to the dot 
capacitively, and thus can modulate the energy of the 
electron localized in the dot with the amplitude eVdot- 
We assume that the leads can be either only dc-biased, 
or an additional ac bias can be applied. 



B. Time-dependent Schrieffer- Wolff transformation 

In the present paper we consider the dot in the Kondo 
regime, U — Ed, Ed ^> Fl,r- Under such conditions, the 
number of electrons on the dot is a well-defined quantity. 
In the limit of infinitely small tunneling, the ground state 
of the system described by the Hamiltonian (Q) is dou- 
bly degenerate due to the spin of the (single) electron 
which occupies the level d. The states with two or no 
electrons on the dot are higher in energy by U — Ed or 
Ed respectively, and are not important for the low-energy 
dynamics of the system. In our paper, we study the ir- 
radiation effects when the applied fields do not drive the 
dot out of the Kondo regime, 



Therefore, the excited states with two or no electrons on 
the dot are to be projected out. This-can be achieved 
by the Schrieffer- Wolff transformation,^ modified to ac- 
count for the time dependence of the parameters of the 
Hamiltonian ([!]). In the present subsection we perform 
this transformation, which finally yields the description 
of the quantum dot system by means of the Kondo Hamil- 
tonian with time-dependent parameters. 

First we move all dependence of the Hamiltonian on 
the applied voltages Vdot, Vdc, U ac to the off-diagonal 
terms. It is achieved by the unitary transformation 



U = exp 



T 1 dt 



E ^(*')4 

k.a.a 



V dot {t')dUa 



(3) 



After this transformation, the Hamiltonian has the form 



<9U. 



H' = Ui/U t -ih— U f 
dt 



E ^l a c kaa + E(-^)4rf CT + uA d A d i 



E 



„Ju + H.c. 



with 



v a (t) = v a exp 



re 



dt' [V a (t')-V Aot (t')] 



(4) 



(5) 



Now we can make the time-dependent Schrieffer- Wolff 
transformation, which is defined by the unitary operator 



W = exp(V) 

with 

v = E {K2(*)(i-^)+4 2 2 (On- 



ce) 



H.c 



} • (7) 



The functions (t) are to be found from the condition 



<9V 

= H v + [V, Hi] - ih— 



(8) 



eV Ac , eUdot, eVac < E d , U - E d 



(2) 



where H v is the part of the Hamiltonian (|j) responsi- 
ble for mixing of electron states in the leads and in the 
dot [the last term of Eq. (Q)], and Hi describes the un- 
coupled dot and leads [the first three terms in Eq. (Q)]. 
The condition (|J) ensures that the resulting Hamiltonian 
WiJ'W 1 ' has no linear -in-v a terms, which account for the 
variations of the number of electrons in the dot. 
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The only difference of the transformation 
the conventional Schrieffer- Wolff transformatior 
time dependence of loj^ . For the static Anderson Hamil- 
tonian, these factors are constant; in our case they are 
functions of time because of the time variations of the 
Hamiltonian (TO). Solving Eq. (§) for toj^ft), we obtain 



-i{Sk-E d )t/h 



(9) 



When the applied ac fields are slow enough, 
hui/Ed, huj'/Ed <C 1, one can solve Eq. (||) in the adi- 
abatic approximation, neglecting the third term in it. 
This approach yields a simplified expression for (t) : 



.(*) 



s d -ey a (t) + ey dot (t) 



(10) 



Here we have also neglected the single-electron energies 
£/c because the Kondo effect is produced by the states 
close to the Fermi level, whose energy is small in com- 

(2) 

parison to Ed- The formulas for wi a (t) are analogous to 
Eqs. (p|) — (|l0[) , only E d must be replaced by E d — U. 

Applying transformation (H)-© '° the Hamiltonian 
([|), we come to the Kondo Hamiltonian 



H K = H a + Hj, 4=E^L 



a kaa. ' 



(lla) 



Hj=Y. (|<W + S lS l aa ,) cl a c k , a , a , , (lib) 

k. a, a 
k , er, a 

where s and S are the spin operators of the electrons in 
the leads and of the electron on the isolated level, respec- 
tively; we assume summation over the repeating indices 
I = x,y, z. In the framework of the Hamiltonian (|lT|), the 
state of the dot is fully described by its spin. The terms of 
the Anderson Hamiltonian ([!]) that are responsible for the 
electron tunneling to/from the dot, and for the Coulomb 
interaction of electrons in the dot have been transformed 
to the term Hj of the Kondo Hamiltonian (fill). This 
term represents exchange interaction between the spin of 
the dot and the electrons in the leads. The Hamiltonian 
(|ll|) operates within the band —Ed<£k < U — Ed, see 
Rcf. |lj. The coupling parameters J aa i (t) are given by 



Jaaf(t)=V a (t) \-W^ a> {t) + W^) al {t) 



C.c. (12) 



The applied bias is accounted for by the time depen- 
dence of Jaa'it) with a a' . The dependence of the 
right-hand side of Eq. ( |l2| ) on the indices k is negligi- 
ble near the Fermi level; therefore we can disregard this 
dependence in the coupling constants J. 

In this paper, we are primarily interested in the irradi- 
ation effects when the applied fields are unable to ionize 
the dot, 



eV dc , eV d ot, eKc < E d , U - E d , 

hoj, hu' <:E dl u - E d . 

Under these conditions, one can use the approximate so- 
lution (|l0|) of the equation (|^), expanding it in powers 
of small parameters eV& t/Ed, eVdot/{U — Ed), etc. For 
simplicity we will consider the cases when the system is 
affected by only one kind of ac field: either ac voltage 
applied to the gate, or the ac bias. 

In the former case, eVdot 7^ 0, eV ac = , we obtain the 
following expression for the coupling parameters: 



Jaa> (t) = Jaa> I 1 + 7 cos ^t] exp 



Vdc,a')t 



where the exchange constants are given by 



J, 



(0) 



irvE d 



E tl 



(U - E d )E d 
U 



(13) 



(14) 



The exponential factor in Eq. jT^ ) is due to the dc bias, 
which produces the phase difference between the elec- 
trons in the left and right lead. The cosine term accounts 
for the applied ac field and stems from the adiabatic vari- 
ation of the electron energy in the dot, Ed + eVdot(i), see 
Eq. (|l(]). The strength of the applied ac field is charac- 
terized by the dimensionless parameter 



2E d - U 
1 = eV **{U-E d )E d ^ 



(15) 



If the ac field is applied to the leads rather than to 
the gate, eVdot = 0, eV ac ^ 0, the expressions for J aa i (t) 
read 



J LR {t) = J^exp 



ieV dc t 



ij' sin(a/i + 4>q) 



(16) 



where J^a 1 IS gi ven by Eq. (|lj). The ac bias creates 
the phase difference between the electrons in the left and 
right lead, and therefore enters the exponent in (fL6|) to- 
gether with the dc bias. The relevant parameter charac- 
terizing the strength of the ac perturbation here is 
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ha' 



(17) 



The variation of the electron energy in the dot with re- 
spect to the leads, E d ± eV r ac (t), see Eq. (|Tc|), generates 
a term smaller by a factor of ~ Tito' / E d , and is neglected 
in Eq. ©. 

In the limit of small amplitude of ac bias, 7' <C 1, the 
expression ( |l6| ) for J aa i (t) may be further simplified by 
dropping the terms of high orders in 7'. Expanding the 
right-hand side of (|l^) in powers of 7' up to the first 
power, we arrive at 
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ieVd c t 



[1 + i-y' sm{uj't + fa)] . (18) 



C. Kondo conductance in equilibrium 



J, 



(0)1 



LR 



J, 



(0) 



RR 



I) 







dh / dt 2 (S j (0)Sk(ti)S l (t 2 ))^ kl 



(22) 



x [t 1 cosieVh/h) + t 2 cos(eVh/h)] 
x G kl (-t 2 )G k2 (t 2 - t^Gk.ih) , 



In the framework of the Kondo Hamiltonian (|ll|)-(|12|), 
two types of tunneling between the left and right leads 
are possible: regular elastic cotunneling [the first term in 
parentheses in Eq. ( |llb| )], and "exchange cotunneling" 
[the second term]. In an act of "exchange cotunneling," 
the simultaneous flip of the spins of the tunneling electron 
and of the dot can occur. In the case of weak coupling 

{v\Jaa'\ u ^ 1)> one ma y a PPly the perturbation theory 
to evaluate the conductance through the dot. It turns out 
that at T — > 0, the higher-order terms of the perturba- 
tion theory series grow, finally making the series diverge, 
signaling the Kondo anomaly. This phenomenon was ex-, 
tensively studied for t he magnetic impurities in metals .u 
In the subsection II C we demonstrate how a similar be- 



havior emerges in the tunneling through a quantum dot. 
The main purpose of this subsection is to present the 
formalism which is suitable for the treatment of a non- 
equilibrium case at hand. For simplicity, we first consider 
the case of no ac field. Effects of the ac field are included 
into consideration in subsequent sections. 

Unlikej-, the conventional treatment of the Kondo 
problem]!! we have to consider the Kondo anomaly di- 
rectly in the conductance, rather than in the scattering 
amplitude. This need emerges from the kinetic nature 
of the problem at 7,7' 7^ 0. To calculate the differential 
dc conductance G(Vd c ), we employ the non-equilibrium 
Keldysh technique in the time representation. In this 
formalism 



G(Vdc 
where 



d 



dV d , 



■<S(-oo, 0)1(0) S(0,-oo)) 0> (19) 



^koL^k'o'R. 



H.c. 



(20) 



is the current operator, and S(i2, t\) is the evolution ma- 
trix determined by Hj. 

In the lowest non-vanishing (second) order of the per- 
turbation theory in the coupling constant J^, , the con- 
ductance of the dot is given by the expression 



G (2) 2 

nh 



-t(0) 
J LR. 



(21) 



The logarithmic divergences appear starting from the 



e^ kl is the antisymmetric unit tensor and Gfc(i) and Gfc(i) 
are the time-ordered and anti-time-ordered Green func- 
tions of free electrons in the leads, given by 



G k {t) 



" /fe)] 
i /(€*)] , 



if 
if 



t > 
t < 



(23) 



with /(£) being the Fermi distribution function. This 
and other terms of the same structure yield the Kondo 
divergence in the conductance. 

If there is no external ac field, the averages 
(Sj(ti)Sk(t2)Si(t3)} are independent of time and equal 

{i/4)ejki- After adding up all the cubic in J^, terms 
in the expression for the conductance G [one of them is 
given by Eq. (p2|)1, summing over the electron states ki, 
and performing the integration over i 2 [see Eq. d22)], we 
arrive at 



>(o) 

J LR. 



dt ■ 



2 irh 
(—t) cos(eVdc</^) 



-7(0) 

J RR 



7(0) 
J LL 



sinh^vrTt/ft) 



ttT 

(T/A,) 2 



-r (24) 



Here 



D a = ^E d {U-E d ) (25) 

is the effective bandwidth0 For the sake of simplicity, 
further we will mostly consider the zero-bias conductance 
Gpoak- In this case, Eq. ( pi| ) yields 



3tt2 ef_ 
2 irh 1 



J, 



(0) 



LR 



J, 



(0) 



RR 



-t(O) 



In 



Do 
T 



(26) 



The results for the finite-bias conductance G^ (Vd c ) with 
eVdc > T can be obtained from Eq. (26) by replacing T 
with eVdc- 

Thus the second [Eq. @] and third [Eq. (|2|)] orders 
of the perturbation theory in the coupling constant 
yield the following expression for the dot conductance 



G 



peak 



4 nh 



4 nh 



7(0) 

J LR 



7(0) 
J LR 



2- ( J R l 



In 



£0 
T 



(27) 



(0) 



terms of the third order in J^ a 
has the following structure: 



A representative term Here we have split the quadratic in contribution 

( fH| ) in two: the one due to the "exchange cotunneling", 
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which entered the first term in Eq. (p7|), and the one 
due to regular cotunneling, which became the last term 
in Eq. @. The cubic in jj®, term in Eq. (§7|) grows 
as the temperature is being lowered, demonstrating the 
Kondo anomaly. The regular cotunneling does not pro- 
duce terms growing at low temperatures and bias, and 
does not contribute to the Kondo effect. Equation ( p7| ) 

is valid while T > T K ~ D Q exp[l/K jf^ + J R °l)]- 

If this condition is not satisfied, then the expansion up 
to the cubic in J^°J, terms is insufficient. At T > Tk, 
the conductance can be derived in the leading logarith- 
mic approximation. The latter consists in the summa- 
tion of the most diverging terms in each order in 3^1' > 
i.e. terms, proportional to [Jlr] 2 [J3 \n(D /T)] n , in 
the series for G. To perform this summation, we modify 
the "poor man's scaling" tcchniquclij In the framework 
of this technique, the electron bandwidth D is gradually 
reduced, and the exchange constants in the Kondo Hamil- 
tonian (|llj) are renormalized to compensate for this band 
reduction, i.e. J^l IS replaced with some J aa > (D) . The 
proper dependence of J aa > on D should be derived from 
the condition of invariance of physical quantities with re- 
spect to the RG transformation. Finally, the renormal- 
ized Hamiltonian with the reduced band width will allow 
for the calculation of the conductance in the second order 
of the perturbation theory in the renormalized exchange 
constants J aa ' ; the resulting expression will be equal to 
the sum of the dominant terms of all orders of the per- 
turbation theory in the initial, bare exchange constants 

j3(d = d ). 

For the non-equilibrium system we consider, the 
Renormalization Group (RG) equations for the exchange 
constants should be derived from the condition of the 
invariance of the linear conductance (or current) under 
the RG transformation, rather than the invariance of the 
scattering amplitudes. In the main logarithmic approx- 
imation which we are going to employ, the (invariant) 
conductance must be evaluated in the two lowest non- 
vanishing orders of the perturbation theory, namely in 
the second and third ones, see Eqs. (^fj), (p4|). The 
Kondo divergence (and, therefore, the renormalization of 
J aa i ) occur due to ex change scattering [the second term 
in braces in Eq. ( lib )] only. Therefore we single out this 
contribution in the term of the second order in Jaa 1 , 



G 



(2) 



3tt 2 e 2 



(28) 



In the third order in the exchange constants, the con- 
ductance is given by Eq. (pi|). The resulting condition 
of invariance of G under the transformation, which cor- 
responds to the "poor man's scaling", has the following 
form: 



dDl exch 
3^ 2 



(D) 



2 nh 



[Jlr] [Jrr + Jll] In 



0. (29) 



Within the accuracy of this equation, when differenti- 
ating the second term, we should neglect any implicit 
dependence on D through the parameters J aa '(D). 

Equation (|2g|), together with Eq. (p8|), yields the equa- 
tion for the evolution of Jlr'- 



AJlr _ Jlr{Jrr + Jll) 
dD ~ V D 



(30) 



The corresponding equations for Jrr and Jll can be 
derived from the condition of invariance of other phys- 
ical quantities under the RG transformation. For this 
purpose, we pick the spin current from the left and right 
lead, 



(S(-oo,0)jW(Q)S(0 > -oo)) 0> 



H 



J, 



k 



(4 



k1a C kfa 



'feia C fe|a 



(31) 
(32) 



which is induced by applying infinitely small magnetic 
field to the leads. The resulting two equations will be 
independent, in contrast to the corresponding equations 
for the charge, because the spin of the dot can vary while 
the charge cannot. 

( s) 

Evaluating 7^ ; in the second and third orders of the 



perturbation theory in J^°J, , similarly to Eq 
differentiating it by D, we arrive at 



<MjiR = ^ Jrr + Jf 
dD 
dJ L L 



LR 



dD 



D 

, ^LL + 3lR, 

D 



), and 

(33) 
(34) 



Equations (|30|), (|33|), and (34) make a complete system, 
which, with the initial conditions 



Jaa'(D ) 

see Eq. (|l|), yields 
Jlr(D) = 



J, 



(0) 



irisE d 



2Vr7T^ 



f 



T L + T R 2uhx(D/T K ) ' 
The Kondo temperature Tr- here is given by 



r R )u 



■ exp 



(35) 



(36) 



(37) 



with 



1L 



L To obtain the pre-exponential factor 



(Tl + Tr)U/it in the equation for Tk one has in fact to 
include the next in J aa '(D) order in the RG equations, 
see Ref. || 

The renormalization should proceed until the band 
width is reduced to T. After that, the current and con- 
ductance can be calculated in the Born approximation 
(p8|) in the renormalized exchange constant Jlr given 



G 



by Eq. (|36|) with D — T. The resulting expression for 
the conductance in the domain T > Tk is 



G 



3^ 



1 



peak 



16 MT/T K )\ 



where 



Gtj = 



(38) 



(39) 



{T L +T R ) 2 

is the conductance of the dot in the unitary limit of tun- 
neling. 

At T 3> Tk, one can expand Eq. into the series 
of powers of J^, \si{Dq/T). The first term of the series 
is the conductance calculated in the Born approximation 
[see Eq. (pq)], the second term yields the lowest order 
Kondo correction given by Eq. (Uq). 

The RG technique can be also used to derive the depen- 
dence of Kondo conductance on the applied dc bias in the 
domain eVdc It Tk, eVd c > T. Starting from Eqs. (p3) 
and (|2|) and proceeding along the lines of Eqs. (p9|)-(p6P, 
we arrive at 

G(V dc ) = ^ 1 Gu . (40) 

16 [\n(eV dc /T K )\ 2 

Thus the Renormalization Group technique (^9|)-(|3^) 
allows one to perform summation of infinite series of the 
perturbation theory in the exchange constants J^^fJ, ■ The 
results obtained in this way are valid in a wider domain 
of parameters as compared to the results of the finite- 
order perturbation theory. The RG technique reveals 
the meaning of the energy scale Tk- The resulting ex- 
pressions (|3g), (fh)|) for physical quantities contain the 
single relevant characteristic of the system, Tk, rather 
than numerous parameters of the Anderson Hamiltonian 
[Eq. ([Taj)] . For example, in Eq. (|3|) the dependence of the 
differential conductance on the applied bias is expressed 
in terms of the dimensionless variable T /Tk- The depen- 
dence of GjGjj on this variable iSpgiven by some univer- 
sal function of any value of T/Tk'p its high-temperature 
asymptote (|8|) is established with the help of RG tech- 
nique. Similarly, the frequency and magnitude of the 
applied ac field may enter into some new universal for- 
mulas for G/Gjj in the form of dimensionless variables, 
being normalized by Tk- The generalization of the RG 
technique which we presented in this section, will allow 
us to check the validity of this conjecture and to establish 
the asymptotes of these new universal dependences. 



III. SPIN DECOHERENCE BY AC GATE 
VOLTAGE 

Now we include into consideration the effects of an 
ac field. As we have shown in our earlier paper Ja the ac 
field can bring decoherence into the dynamics of the dot's 
spin, thus affecting the Kondo conductance. We start our 
study of the irradiation-modified Kondo anomaly from 
the consideration of the decoherence. 



A. Mechanisms of spin decoherence 



In terms of the Anderson Hamiltonian (|la| ) the loss of 
coherence by the dot spin occurs when an electron leaves 
the dot and another electron, with the opposite spin, en- 
ters it. If the frequency of the applied ac field is large 
enough, Hto > Ed, U — Ed, this process can consist of two 
real processes: the dot gets ionized by the ac field, and 
then an electron from a lead enters the dot to fill the va- 
cancy. Alternatively, an extra electron can be put in the 
dot, and then the electron which was initially present in 
the dot leaves it. 

In the present paper we deal with a more subtle case, 
when the applied ac field is unable to ionize the dot. In 
this case the dot can still change its spin, even at zero 
bias, by means of the "spin-flip cotunneling," which is 
shown schematically on Fig. [l]. In the course of this pro- 
cess, an electron, which interacts with the dot spin [see 
Eq. (p"l|)], absorbs a photon and hops to a state above 
the Fermi level, while the spin of the dot flips. In terms 
in the Anderson Hamiltonian ([la]), this process cannot 
be described as two separate real processes. Instead, the 
change of the dot spin occurs as a single process when a 
state with two or no electrons in the dot appears only as 
a virtual intermediate state. 



(0 



(f) 



~ha> 



FIG. 1. Spin- flip cotunneling: absorbing a photon, an 
electron hops from a state below the Fermi level to a state 
above the Fermi level; the spins of the electron and of the dot 
flip due to exchange interaction between them. 

The rate of the spin-flip cotunneling can be calcu- 
lated with the help of the Kondo Hamiltonian given by 
Eqs. (H), and (0). In the case of weak modula- 

tion, 7 <C 1, see Eq. (|16[), it is sufficient to account for the 
single-photon processes only, and use the reduced form 
of the Hamiltonian, given by Eqs. (11), (p~3|) - (p~5|) - An 
infinitely small dc bias, needed for measurements of the 
linear conductance, does not affect the rate of spin-flip 
cotunneling. Therefore in this subsection we set Vd c = 
for the sake of simplicity. 

Applying the Fermi Golden Rule we obtain 



— = — TVjJ 
t 8tt 



E d 



(41) 



7 



where 7 is given by Eq. (|l5|). 

The amplitude of inelastic transitions yielding Eq. (El]) 
was evaluated in the lowest-order perturbation theory. 
It corresponds to the first order in the amplitude of 
the ac perturbation, and the zeroth order in the time- 
independent (at Vdc = 0) part of the exchange interaction 
( |llb| ). Accounting for the terms of higher-order in this 
time-independent part renormalizes the amplitude of the 
inelastic transition (which is still linear in the amplitude 
of the ac field). Similarly to the calculation of the conduc- 
tance, we intend to collect the leading logarithmic terms 
in the renormalization of the inelastic transition ampli- 
tude. This can be accomplished by the renormaliz ation 
group (RG) transformation, described in subsection [IC . 
The transformation reduces the electron band width D 
and renormalizes the matrix elements J aa ' of the Kondo 
Hamiltonian ([ll]) to account for this band reduction. Fi- 
nally one can calculate the decoherence rate in the second 
order perturbation theory in renormalizcd J aa i ; the re- 
sult given by such a treatment equals the sum of infinite 
perturbation theory series in the initial Hamiltonian. 

The RG transformation starts from the bandwidth 
D = Do given by Eq. (p5|), and the initial matrix ele- 
ments 



J'aa' I D=ZJo = ^ [1+7 COS Wt] 



(42) 



cf. Eq. (|l3|). While the width D of the band exceeds 
hcu, the time dependence of the Hamiltonian matrix ele- 
ments ( |42] ) can be treated adiabatically, i.e. time t in the 
right-hand side of Eq. ( f42| ) can be considered as just a pa- 
rameter. The RG equations, derived from the condition 
of invariance of physical quantities under the transfor- 
mation, have the same form (|30|), (Q), and (|34|). The 
transformation must be stopped when the bandwidth is 
reduced to the values of the order of the frequency Hlo of 
the applied ac field. Expanding the solution of the RG 
equations (|30|), (|33|), and j3^ ) with the initial condition 
(H2) in powers of 7 up to the first power, we obtain 



2-\/r Q r Q / 



1 



nE d 1 



1 + 7 



r L + T R \n{huj/T K ) 



cos Lot 



T K (t) = n\l exp 



(Tl + T R ) 



with 



Ed(t) = E d (l +7coswt) 



(46) 



(47) 



(43) 



The Fermi Golden rule applied to the Hamiltonian (ft 
with J aa i given by Eq. ( |43| ) yields the following expres- 
sion for the decoherence rate: 



37T TlLJ 



1 



ST, 



K 



cf. Eqs. and (§J). 

One can see that the amplitude of the ac field enters 
Eq. (0) through the dimensionless parameter 5Tk/Tk- 
The value of this parameter, in principle, can be directly 
measured. Representation of Ti/tTk in terms of STk/Tk 
will allow us to build a universal description of the ac field 
effect on the Kondo conductance. 

As we mentioned before, the spin-flip cotunneling is 
essentially different from the dot ionization with its sub- 
sequent filling. During a process of spin-flip cotun- 
neling, the ionized dot exists only as a virtual state. 
Therefore, the spin-flip cotunneling persists at frequen- 
cies minlii^, U — E d } > hcu, leading to the decoherence 
of the dot spin state without ionization. 



Effects of spin decoherence on Kondo 
conductance 



As we have just shown, the external ac field is able 
to flip the dot's spin. Therefore, in the presence of the 
ac field, the averages of the type (Sj(ti)Sk(t2)Si(t 3 )) no 
longer equal (5^(0)^(0)^(0)) = TtSjS k Si = (i/fyjkl- 
In the limiting case \t m — t n \ 3> t (m 7^ n), the orienta- 
tions of the dot spin at t = ti, t%, £3 are independent of 
each other, because of the spin-flip cotunneling, and one 
has 

{SjWSkfoWfa)) = (Sj(ti)){S k (Jt2)){Si(t 3 )) = . 

At finite time intervals \t m — t n \, the spin correlator de- 
cays exponentially, with the spin-flip cotunneling rate be- 
ing the characteristic decay rate: 

03{ti)S k (t2)Si(ti)) = {i/±)e jkl exp{-t max /T) , (48) 
t max = max{|ii - t 2 \, \t 2 - t 3 \, \ti . 

Equation (|f8|) can be derived using the formalism of 
equations of motion. In the framework of this formal- 
ism, Eq. (f48|) appears as the solution of the equation 



tT k 32 T K [\n{huj/T K )] A 
Here we introduced the relative amplitude 
5T K irE d 



T, 



K 



7i 



r, 



(44) 



(45) 



_d_ 



(t 3 ,h) (ifi^S,]) S(ti,*a)-§fcS(ta,*s)^>o (49) 



of adiabatic variations of the "time-dependent Kondo 
temperature". The latter is defined by 



= (S 



where S(i, t') is the evolution matrix determined by Hj. 

Expanding S(t,t') in powers of ^jf^, up to the first 
power, we arrive at 



8 



■^(Sj{h)S k (t2)Si(t 3 )) 
= - [6(t 3 - h)6(t 2 - h) - 0(h - ta)6(ti - h)} 



(50) 



where r is given by Eq. ( |4l] ) . Equation ( |50"| ) with r given 
by Eq. (Q) can be obtained by expanding the evolution 

matrix S(t, t') up to the second pow er in y j^, , and using 



the RG technique described in Sec. II C 



The leading effect of the irradiation is in cutting off 
the logarithmic divergences in the time integrals of type 
(E3) . One can easily see that with the time-decaying spin 



correlation function 



G 



(3) 
peak 



is finite even at T —> 0: 



2 e 2 



3tt 



2 TT% 



7(0) 



J, 



(0) 



RR 



-r(O) 



hi 



or 



(51) 



As we have shown in Sec. Ill A, the spin decoherence by 



external irradiation does not require ionization of the im- 
purity level, and therefore exists at frequencies below Ed, 
U-E d . 

The effect of the irradiation on the Kondo conductance 
is not analytic in the intensity of the ac field. It cannot 
be obtained by a finite-order perturbation theory in 7 
in the formula (|l9|). To obtain Eq. ( filf ) directly from 
Eq. (^) using the perturbation theory series in jj^, , 
one would need to add up all the terms proportional to 

The finite-order perturbation theory [Eqs. ( |2l|) and 
(pl|)] can be used to evaluate the Kondo conductance only 
if the decoherence rate %/t is much larger than the Kondo 
temperature Tk- At lower decoherence rates we have to 
take into account terms of all orders in J^a>- ^ can be 
done by means of the Renormalization Group technique 
described in Sec. 



IIC 



One RG equation is to be de- 
rived from the condition of invariance of the conductance, 
given by the second and third orders of the perturbation 
theory in J^, [Eqs. (§|) and (|l|)], similarly to Eq. ©. 
The other two RG equations can be obtained using the 
requirement of invariance of the spin current @ under 
the RG transformation. The resulting set of equations 
coincides with the one given by Eqs. (pp|), (|33|), and (34). 
When the decoherence rate exceeds the temperature T, 
the RG transformation must be stopped when the band- 
width D reaches %/t rather than T. Then the linear 
conductance can be evaluated in the second order per- 
turbation theory in the renormalized exchange constants 
J LR: given by Eq. © with D = %/t: 



G 



3tt 2 



1 



peak 



16 [\n(%/TT K )Y 



:G 



u ■ 



(52) 



Here %/tTk is determined by Eq. ([44]). Equation ( |52; ) 
is the central formula of this section. It defines the con- 
ductance of the quantum dot as a function of two di- 
mensionless parameters: %uj/Tk and STk/Tk [Eq. (Eq)]. 



The region of validity of Eq. (|52j) is determined by the 
condition 



tT, 



> 1 



A" 



(53) 



and corresponds to the regime of strong suppression of 
the Kondo effect by the external ac field. At fixed 
strength of the ac field the spin-flip rate (|4l]) decreases 
with the decrease of ac field frequency uj. Correspond- 
ingly, the peak conductance ( |5^ ) grows. The crossover 
from weak to strong [G ~ Gjj] Kondo effect occurs when 
%/t ~ T K . Equations ([id]) and (Q show that this value 
of %/t is reached while %lu/Tk S> 1. 



IV. WEAK SPIN DECOHERENCE 

In this section we consider the regime of "intermedi- 
ate suppression" of the Kondo effect by the ac radiation. 
By "intermediate" we mean that the decoherence is rel- 
atively weak 



%/t < T K , 



(54a) 



and the formula ( p2| ) is no longer valid, but the frequency 
is still sufficiently high 



Tk < %u> 



(54b) 



so that the RG result (44) for the decoherence rate holds. 



In this regime, the formation of the many-body state is 
not suppressed, because of Eq. (54a). However, Eq. ( |54b| ) 
allows for sudden spin flips. The complicated nature of 
the many-body state hampers the quantitative consider- 
ation of this regime, and we limit ourselves to qualitative 
analysis. 

When the many-body Kondo resonance is fully formed, 
the conductance of the dot equals G\j [Eq. (gjj)] and cor- 
responds to the unitary limit of tunneling through the 
dot. An act of spin flip destroys the many-body state, 
and the conductance drops substantially below the value 
given by Eq. (39). The time necessary for the mauy-body 
state to be restored equals approximately %/Tk Ji3 There- 
fore, the fraction of time which the system spends in the 
highly-conducting (G w Gjj) state equals approximately 
1 — clH/tTk, where a ~ 1. The resulting time-averaged 
conductance of the dot can be estimated as 



G 



peak 



1 -a- 



tT k 



G, 



(55) 



The rate %/t of the spin- flip processes here is given by 
Eq. ([iil). Under conditions (|54|), parameter a does not 
depend on the characteristics of the ac field. The value 
of a should be found from the quantum-mechanical prob- 
lem of evolution, which starts with a state "prepared" by 
the flip of the impurity spin, and results eventually in the 
re-formation of a Kondo polaron. Our qualitative treat- 
ment of the regime (B4J) does not allow us to find the 
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exact value of the universal coefficient a, which however 
can be found from a numeric calculation. 

At the upper limit of applicability, %/t ~ Tk, the peak 
conductance given by Eq. (^) matches the result (p^). 



V. LOW-FREQUENCY AC FIELD: ADIABATIC 
APPROXIMATION 



In Sections [I C and III we considered the case when 



the isolated spin is only weakly screened by the many- 
electron state formed around it. The complete screen- 
ing was suppressed either by relatively high temperature, 
T > Tk, or by large bias eVdc > Tk, or by the decoher- 
ence. In the case of Sec. [V, the spin-screening cloud is 
able to form; however, the spin flips, produced by the 
irradiation, occasionally destroy this many-body state, 
thus reducing the conductance. 

In this section, we consider the case of low frequencies 
of the ac field hu <C Tk, when the energy of a photon 
is insufficient to flip the dot's spin in the fully formed 
many-body Kondo state. For the irradiation to be the 
leading cause of the deviation of the conductance from 
the unitary limit, we suppose the temperature and bias 
to be also low, T, eVdc <C Tk- The RG technique we 
used before is not applicable in this regime. Therefore 
we need another approach to evaluate the conductance 
of the quantum dot system and the effects of the external 
irradiation on it. 

The required approach is provided by the scaling the- 
ory of Nozieres and Blandin.EZl This theory states that 
the renormalization-group transformation, whose initial 
stage was described in Sec. 



[I C , can be continued, and 



finally leads to a fixed point. At the fixed point, the sys- 
tem exhibits Fermi-liquid behavior, and its Hamiltonian 
has a relatively simple form.tS This fixed-point Hamilto- 
nian can be used to study the prope rties .of the Kondo 
system at low temperatures, T <C Tk-.Oi13 Mapping the 
quantum dot system in the Kondo regime onto the regu- 
lar one-channel Kondo problem, we can employ the fixed- 
point Hamiltonian to evaluate the dc current through the 
dot induced by the applied bias. 

The external ac field disturbs the many-particle state 
formed near the isolated spin, leading to the deviations 
of the system behavior from that dictated by the (static) 
fixed-point Hamiltonian. In this section we study the 
case when the frequency of the field is low (hu> <C Tk), 
so that the many-body state is not destroyed but rather 
adiabatically varied by the ac field, as the level in the 
dot goes up and down [see Eq. dig)]. Then the cur- 
rent through the dot can be evaluated with the help of 
the fixed-point Hamiltonian with time-dependent param- 
eters. 

Now we map the problem of transport through the dot 
onto the regular scattering problem. For this purpose, it 
is convenient to use the basis of s and p scattering states 
rather than that of the left-lead and right-lead states. 
These two bases are connected by 



4o ) = £> c koL + r\c kaR 
vl 



>) 

her 



where £ 



-VCkcrL + £c kc rR , (56) 
vr 



The p-states are decoupled from the dot, so the dot-lead 
coupling term in the Anderson Hamiltonian ( |la| ) has the 
form 



The initial basis Cka-a is composed of the states residing 
entirely in the left or right lead, which is convenient for 
the problem of two leads connected by a weak link, when 
the inter-lead tunneling is to be considered as a pertur- 
bation. In terms of incident and reflected/transmitted 
waves, these states correspond to the waves incident from 
one of the leads to the dot and completely reflected back 
to the same lead. Therefore the s-waves of Eq. d56j), 
which enter the new basis, have the scattering phase 
equal to tt/2. 

Making the Schrieffer- Wolff transformation, we arrive 
to the regular Kondo problem, which at low tempera- 
tures can bej-studied with the help of the fixed-point 
Hamiltonian.lia Under these conditions, the s-wave elec- 
trons, interacting via the isolated spin, form the screening 
cloud. This many-body state still has the Fermi-liquid 
properties, though its scattering characteristics are dif- 
ferent from those of just an isolated spin. One of the 
principal differences is the shift of the scattering phase 
by tt/2 for the states at the Fermi level.Eil This suggests 
another change of the basis for the sake of convenience: 
from s-waves having the scattering phase equal to 7r/2, 

a icr' t° those with the scattering phase tt. 

(s) 

The formal relation between the a k J and the new basis, 
which we denote b ka , is given by 



dxe lkx ^ a {x) , 



in / dx' g(x') 



^cr{x) = exp 



(57) 



where g{x) is an arbitrary function obeying 
Jf° dx' g(x') = 1. The "coordinate" x was introduced 
for convenience to separate the incoming and outgoing 
parts of the scattering states, which correspond to nega- 
tive and positive values of x respectively. 

Before the scattering region [x — > — oo in Eq. ([571)1, 
the wave functions of the states b ka and ail coincide. 
Therefore the states 

J in ) — Ch m „W 



„( in ) 



7R = vb k ,+^ (58) 



represent waves incident from a left or right lead only. 
Passing the scattering region, the wave function of the 
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state bka acquires an extra phase of 7r as compared to 

(s) 

that of a k( ^ . Then one can see that the states 

(out) (h (p) (out) _ , e (p) ,r Q \ 

c kaL ~ -?°fc CT - V a ka , c kaR - -i]b ka + £a krJ (59) 

have an outgoing wave only in the left or right lead re- 
spectively. The current operator in terms of these states 
equals simply 



I(V) 



E/„(° ut )t„(out) 
y-kaL u kaL 

k.a 



(out)f (out) 
^koR kaR 



(60) 



(s) 

The fixed-point Hamiltonian in the basis bka, has 
a relatively simple form:E3 



k<7 



ka 



V F 



VT: 



K 



fci feacr 



kxkikzki 



where : ... : denotes normal ordering. The spectrum 
of electrons is linearized, = v F k, since the reduced 
bandwidth is of the order of Tr- <C e F ; The Kondo tem- 
perature Tk is the only energy scale of the fixed-point 
Hamiltonian (|6l|). 

The third term in Eq. (pi]) determines the phase shift 
which a quasiparticle acquires as it passes through the 
dot. This shift is energy-dependent: it equals 7r at the 
Fermi level, pas we discussed before; and 7r + /Tk in 
general caseJii In terms of waves incident from the left or 
right lead, such behavior of the phase shift is analogous 
to those in tunneling through a resonant state tied to 
the Fermi level. The fourth term in the Hamiltonian de- 
scribes the interaction of the quasiparticles of the Fermi 
liquid at the fixed point. The p- waves are not affected by 
the Kondo screening, so the Hamiltonian for them has 
the same form as the one given by Eqs. (|la|), (|56|). 

Using the Hamiltonian (|6~lj), we can rewrite the current 
operator (60) in the form more convenient for the further 
calculations: 



f P , 



E (out)t (out) 
^kaR c kaR 

k.a 



(62) 



The first term in Eq. (|62| ) is the current that would flow if 
all the states were able for resonant tunneling through the 
dot; the scattering between the left- and right-incident 
species (which is just backscattering when £ = r)) re- 
duces the magnitude of the current, and is accounted for 
by the second term. 

To evaluate the conductance of the dot, we employ the 
Keldysh technique (fi9|), treating the last two terms of 
Hamiltonian (pl|) as a perturbation. At infinitely small 



temperature and bias, the current through the dot is 
transferred by electrons at the Fermi level. The trans- 
mission coefficient for these electrons equals (2£?7) 2 , i.e. 
the second (backscattering) term in the current opera- 
tor ( |62| ) yields zero. Therefore the dot at this conditions 
has maximum conductance, G = Gu = (e 2 /Trh)(2rj£ i ) 2 . 
At finite temperatures, the electrons which facilitate the 
current are spread within a strip of width T near the 
Fermi level. The departure of the particle energy from 
the Fermi level in the system (^Tj) leads to the deviation 
of its scattering phase from tt, i.e. from the resonance. 
Therefore the conductance in this case will be lower than 
Gu. Indeed, substituting Eqs. @ and @) into Eq. @ 
and employing the second order of the perturbation the- 
ory in the last two terms of Hamiltonian (^), we arrive 
at 



K 



v F v 3 T% ^ 
K feife 2 fc 3 



df(v F ki) 



dki 



df(v F k) 
dk 

f(v F k 2 ) 



x [l-f{v F k 3 )]f[v F (ki-k2 + k 3 ) 

2" 



1-7T* 



T 



G 



u 



(63) 



where /(e) = l/[exp(e/T) + 1] is the Fermi distribution 
function. One can see from Eq. ([33]) that the conductance 
of the quantum dot system at low temperatures decreases 
with growth of the tempjecature. This behavior has been 
observed experimentalljO'El and is analogous to the de- 
crease of the resistivity in a regular Kondo system (bulk 
metal with magnetic impurities). 

The differential conductance of the dot at finite bias 
Vdc, with T -c eVdc *C Tk, can be derived analogously 
to Eq. (|63|). The resulting formula 



G(V dc ) 



3 / eVdc 



T 



K 



Gr 



(64) 



shows that G(Vd c ) decreases with the growth of the bias 
applied to the dot. 

Slow (Tilj <C Tk) ac field results in adiabatic time- 
dependence of Kondo temperature, see Eq. (E&). The 
time-dependent part of the Hamiltonian (pH ) with 
1/Tff(i) given by Eq. ( |46| ) accounts for the interaction 
of quasiparticles with ac field. To consider this part of 
Hamiltonian in the conventional terms of electron-photon 
interaction, we expand 1/Tr-(£) in Fourier series: 



1 — \^ \ c mu>t 



T K (t) 



n *K 



(65) 



After an act of photon absorption, a quasiparticle facil- 
itating the current is transferred from the Fermi level, 
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i.e. away from the resonance. As a result, at low tem- 
peratures the ac field must reduce the conductance of a 
quantum dot in the Kondo regime. At Uu> <C Tk, the 
conductance can be calculated in the second order of the 
perturbation theory in the time-dependent part of the 
Hamiltonian. Substituting Eq. (|65| ) into Eq. (|6l|), and 
then using the Keldysh formalism ([l9]) to evaluate the 
conductance we arrive at 



break down the universal description of the problem, and 
the Kondo temperature remains the only relevant energy 
scale. We have shown that the effect of the irradiation is 
described by two dimensionless parameters Tioj/Tk and 
STk/Tk oc Vdotj where 6Tk has the meaning of the adi- 
abatic variation of the Kondo temperature under the in- 
fluence of ac modulation, see Eq. (fl5|). 

At sufficiently large frequencies to of the ac field, when 
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where for simplicity we set temperature to zero. Trans- 
forming Eq. @ back from l/T^ to l/T K (t) [Eq. @], 
we finally obtain 
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Gu 



where ... denotes averaging over the period of variation 
of T K (t), and ST K /T K is defined by Eq. @. 

T he sin gle-photon decoherence processes described in 
Sec. [II A do not occur in this regime, because the energy 
necessary to flip the dot's spin is increased by its interac- 
tion with the screening "spin cloud" in the leads, and is 
of the order of Tk ~> Tuu. The rate of the spin flip due to 
many-photon processes is exponentially small in Tk/Tilu. 



VI. SCALING FORMULA FOR THE 
CONDUCTANCE 



Tito 32 [\u(ST k /Tk)} 
Tk~ > 3^ [STk/Tk] 2 



(68) 



the rate Tl/t of the spin- flip cotunncling exceeds the 
Kondo temperature Tk ■ The spin- flip cotunneling brings 
decoherence into the spin dynamics of the dot, destroy- 
ing the Kondo resonance. Small lifetime of the Kondo 
resonance leads to a significant suppression of the Kondo 
effect, see Sec. [II B. The dependence of the zero-bias dc 



conductance G pea k of the dot on the power and frequency 
of the ac field is given by Eqs. (^2|) and (fl4|). 

Upon lowering the frequency lu, condition ( |68[ ) breaks 
down, and Ti/t becomes smaller than the Kondo tem- 
perature. Under such conditions, the strong suppression 
of the Kondo conductance is not possible. However, the 
conductance still may deviate from the unitary limit Gu- 
The violation of the condition ( |68| ) occurs while Tilu still 
exceeds Tk- The zero-bias conductance in this regime 
can be estimated by Eq. (|55|) and (p4|). 

At frequencies below the Kondo temperature, the ac 
field is unable to flip the spin of the dot, and the spin- 
flip cotunneling does not occur. In this regime, the ac- 
driven deviation from the unitary limit is small and can 
be accounted, within the framework of the Fermi-liquid 
description.^ The main role of the ac field is to scatter 
the conduction electrons, transferring them to the ener- 
gies off the Fermi level. These scattered electrons miss 
the Kondo resonance, which is tied to the Fermi level. It 
produces the small deviation of the dc conductance G pea k 
from the unitary limit, see Sec. 0, Eq. ([37]). 

The results obtained for these three regimes match 
each other on the corresponding limits of applicability. 
It allows us to piece together the dependence of G pea k on 
STk/Tk and Tilu in a broad frequency range, see Fig. ^|. 



In this section we summarize the results obtained in 
Sees. pTj|-|v| for the effect of the periodic modulation of 
the dot's potential on the Kondo conductance. 

In the absence of ac irradiation, the quantum dot sys- 
tem is described by a number of physical parameters, 
see Eqs. (|TT|) , (|I4|). However, in the Kondo regime all 
these parameters combine into a single relevant energy 
scale, Tk, see Eq. (|37j), controlling the behavior of the 
system, see, e.g., Eqs. (p8|), and (p3[). The periodic mod- 
ulation Vdot cos tot of the dot potential adds two more 
parameters to the initial Hamiltonian (|i~i|), and, most 
importantly, drives the system into a non-equilibrium 
state. Surprisingly, such a drastic perturbation does not 
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FIG. 2. The zero-bias Kondo conductance G pea k of a 
quantum dot monotonically decreases with the increase of 
the frequency hui of modulation of the dot potential. The 
plot shows the dependence of G pea k on hui at fixed amplitude 
eVdot oc 5Tk/Tk of the modulation. The conductance in the 
unitary limit Grj is given by Eq. (^9|). The value of Gpc^/Gu 
at Ti/t = Tk is denoted by A < 1; this value can be found 
from numeric calculation. 

This dependence allows us to conjecture that at small 
STk/Tk the conductance can be cast in the following 
form: 



G 



peak 



Gu 



F 



ST, 



K 



)"'( 



(69) 



with two universal functions F(x) and f(y). Each of the 
functions depends on only one variable; they have the 
following asymptotes 
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(70) 



(71) 
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The numerical parameter a ~ 1 is introduced and dis- 
cussed in Sec. IV. When fwj ^> Tk, the argument of 



function F has the meaning of dimensionless decoherence 
rate U/tTk- 



VII. 



CONDUCTANCE SUPPRESSION BY AC 
BIAS 



In sections [II -VI we considered the effects of modu- 



lation of the dot potential on the Kondo conductance. 
In the present section we study the conductance in the 
system where the ac field is applied to the leads, thus cre- 
ating alternating bias V ac . The parameters characterizing 
such a held are the dimensionless amplitude eV ac /TK and 
frequency Uui'/Tk- 

First we consider the case of "fast" ac bias, %u>' T^> 
maxjTft-, eV^c}. Under these conditions, the ac bias af- 
fects the Kondo conductance through the decoherence of 
the dot's spin. The dependence of the corresponding de- 
coherence rate %/t' on the amplitude and frequency of 
the ac bias can be calculated with the help of the Renor- 
malization Group technique which we used in Sec. II C , 
III A, and IIIB. The resulting expression reads 
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(72) 



Note that, i n cont rast to the ac modulation of the gate 
voltage (Sec. Ill A ), in the case of ac bias the rate of deco- 
herence decreases with the growth of the field frequency 



The parameter 
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e 2 /irh (T L 



characterizes the asymmetry in the dot, and emerges in 
the expressions for quantities associated with electron 
transfer between the leads. 

When Tl/t'Tk > 1, the conductance can be eval uated 
by means of the perturbation theory, see Sec. |II Cj . The 
decaying function (Sj(ti)Sk(t2)Si(t 3 )), which enters the 
terms of the perturbation theory, provides the large-time 
cut-off for the integrals in equations of the type (B2) . The 
derivation of the expressions for t he conductance is iden- 
tical to the one given in Sec. IIIB, cf. Eq. (52). The final 
formula reads 
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(73) 



with % I r'T K given by Eq. @. 

At smaller amplitudes, Tl/t'Tk < 1, the ac bias is 
unable to suppress the formation of the Kondo many- 
electron state. For this case we may repeat the reasoning 
of Sec. IV. As the result, we obtain 
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(74) 



i.e. the Kondo conductance is only weakly suppressed. 

In the opposite limit of slow variations of bias, Tiuj' <C 
max{TR-, eV^c}, one can use the adiabatic approximation, 



Gpeak = G(Va C COS Uj't) 



(75) 



Here G(V) is the differential dc conductance at finite bias 
V, and 777 denotes averaging over the period of variation 
of ac bias. For eV/Tjc -C 1, the conductance is given by 
Eq. (0). Substituting it into Eq. (|75|), we obtain 
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In the opposite case, eV&c/TK 3> 1, we obtain using 
Eqs. © and @: 
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Figure [| shows the possible regimes of the ac bias effect 
on the Kondo conductance of a dot. At 



Tuv <§C maxjeVac, Tk} , 



(78a) 



the peak conductance depends only on eV ac /TK, see 
Eqs. (|76|), <^n\). In the opposite case of high frequencies, 
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(78b) 



VIII. SATELLITE CONDUCTANCE PEAKS 



the peak conductance depends only on T i/t'Tk, see 

Eqs. @, 0. Thus in both regions ( |78a| ) and ( |78b| ), 
Gp C ak is a function of a single variable. However, the 
corresponding variables arc different in the two regions. 
Therefore, at the crossover between these two frequency 
domains, the peak conductance G vea k(Tiw' /Tk, cV^c/Tk) 
can not be cast into a simple form of a single-variable 
function. 






FIG. 3. The regimes of the ac bias effect on the Kondo 
conductance. The solid line is the border between the two do- 



mains (78a) and (78b). The dashed line separates the regimes 
of weak (below the line) and strong (above the line) suppres- 
sion of the Kondo conductance in each of these two domains. 

It is instructive to consider the peak conductance as 
a function of the frequency u>' at a hxed field amplitude 
eV ac . At small amplitudes, eV ac /IV <C 1, the suppres- 
sion of the Kondo effect is weak at any frequency. At 
stronger fields, eV ac /Tx 3> 1, the Kondo effect is sup- 
pressed far below the unitary limit at low frequencies. 
The height of the zero-bias peak grows with the increase 
of the field frequency u>', see Fig. |L The description of 



the crossover between the regimes (78a) and (78b) is de- 
veloped in Appendix [a|. 
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FIG. 4. The suppression of the Kondo effect with ac bias 
weakens, as the frequency hco' of the bias grows. The plot 
shows the dependence of the dc zero-bias conductance G pea k 
on hoj' at fixed bias amplitude V ac ^> Tk /e. 



In Sections III -VII of this paper we mostly concen- 
trated on the effects of the ac field on the zero-bias Kondo 
conductance G pca k. In this section we study how the ac 
field modifies the finite-bias differential Kondo conduc- 
tance G(Vd c ). 

Without an ac field, the dependence of the differen- 
tial Kondo conductance on Vd c is given by Eq. (64) for 
eVdc < T K and by Eq. (jit]) for eV dc > T K - As we have 
seen from the previous sections, the ac field reduces the 
height of the zero-bias peak, G pca k = G(Vd c = 0). 

Another effect of external irradiation on the differen- 
tial conductance G(Vd c ) is in producing satellite peaks at 
eVdc = infku. If an external ac field is applied, then, at 
eVdc — ±nfko, a tunneling electron can hop from a state 
at the Fermi level in one lead to a state at the Fermi 
level in the other lead, emitting or absorbing n photons. 
Thus at finite bias the external irradiation can effectively 
put a tunneling electron into zero-bias conditions, and 
the Kondo anomaly in the conductance is revived. The 
heig ht of these peaks can be calculated from the formula 
dl9| ) similarly to Eq. (|24|). At low enough irradiation 
level, 7 <C 1, it is sufficient to consider only one-photon 
processes, accounted for by the Hamiltonian (pi]) with 
the coupling constants J aa '{t) given by Eqs. (|l3[)— (p~5|) . 
In this approximation, we will be able to describe the 
first pair of satellite peaks, which emerge next to the 
main, zero-bias, peak. The resulting correction to the 
conductance at e|Vd c | close to Tiw has the form 
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(79) 



When eV 7^ zLTluj, the cosine function cuts off the 
logarithmic divergence here. However, when eVd c — > 
± hui, the cosine factor becomes essentially constant [cf. 
Eq. ( pi| ) at Vdc — ► 0], and the differential conductance 
has a peak again. At T — > 0, the height of the satellite 
conductance peak is determined by the spin decoherence 
rate h/T s th We must mention that r st i may be signif- 
icantly shorter than r given by Eq. (|4l|). The time r 
characterizes the spin decoherence at zero bias, whereas 
the satellite corresponds to a finite bias eVdc = ±fru>. In 
the latter case, the spin decoherence occurs mostly due 
to the tunneling of electrons through the dot (see Fig. ||, 
and also Ref. O). The rate of this process at Tiu 3> Tk 
is given by 



T st i 2tt ^ £2 
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To extend the result ( |80| ) to lower frequencies hui > Tk, 
we employ the RG technique. As the result, we obtain 
the decoherence rate as a function of the universal pa- 
rameter Tius/Tk'- 
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FIG. 5. At finite bias, the coherence of the spin state of 
the dot can be lost in an act of cotunneling, when an electron 
leaves the dot to the lower voltage lead, and another electron 
enters it from the higher voltage lead. 

Using the RG technique, we derive the formula for 
G s ti.±(Vdc), which is the contribution of one-photon pro- 
cesses to the differential dc conductance: 
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(82) 



One can see from Eqs. (|8l|)-(|82|) that G B t\ depends on 
the parameters of irradiation only through the universal 
variables STk/Tk and Huj/Tk- The details of derivation 
of Eq. (|8^) are given in Appendix |b[ 
The full expression for G(Vd c ) reads 

G(Vdc) = G ma i n (Vdc) + G M ,+ (V dc ) + G st i,_(Vd c ) . (83) 

Here G ma i n accounts for the tunneling through the dot 
without absorption/emission of photons and is responsi- 
ble for the zero-bias Kondo peak. The other two terms 
in Eq. ( |33"| ) describe the satellite peaks in G(V dc ). 

As the criterion for resolution of the satellites, we 
adopt the requirement that the function G(V dc ) must be 
non-monotonic on the sides of the zero-bias peak. This 
requirement can be reformulated as: 



ST 



K 



Tj 



> 1 , 



K 



> 1 



K 



(84) 



ij) we used Eq. (|8§) 
Such a form of the 
should be used be- 



In the derivation of the condition ( 
with G main (Vdc) given by Eq. @ 
elastic Kondo conductance G ma i n 
cause at e|Vd c | ~ S-w it is suppressed mainly due to 
the finite bias rather than due to the decoherence, since 
foj) » 7tyV s ti » U/t. 



The above consideration was performed for the ac field 
applied to the gate. The case of ac bias is can be con- 
sidered similarly. The third-order perturbation theory 
result for the shape of satellite peak may be obtained 
from Eqs. @ and @ by replacing 7 2 with (-/') 2 . The 
RG treatment yields 
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The condition for the satellite peaks to be clearly visible 
takes the form 



hoj' 
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Conditions ( |84J) and (p6|), together with Eq. 
demonstrate that upon the increase of the amplitude of 
the ac field, the zero-bias peak is suppressed first, and 
only after that the satellite peaks may become distin- 
guishable from the background conductance. 



IX. COMPARISON WITH EXPERIMENT 

In this section, we discuss possible ways of comparison 
of the theory presented here with experiments. Ideally, 
one should measure the dc conductance through a dot 
while applying the ac bias to the gates or leads in a con- 
trollable way, and varying its frequency u> and amplitude 
in a broad range. Our results predict that the data ob- 
tained at various values of Tk should be scalable, when 
using the proper dimensionless variables, Eqs. (|70|), ([71]), 
©, dfl), ©, and @. We predict also that at a fixed 
magnitude of ac field, the suppression of the Kondo ef- 
fect becomes more severe with the increase of u>, if the 
applied field modulates the gate potential; the depen- 
dence on uj in the case of ac field applied to the leads is 
opposite. If lu significantly exceeds Tk, see Eqs. ( p4] ) and 
(jSq), observation of "satellites" at eV — hto of the main 
Kondo singularity may become possible. The appear- 
ance of even small side peaks though is associated with 
a strong suppression of the zero-bias Kondo singularity. 

Presently, there is only one experiment aimed at obser- 
vation of effects of irpadiation on the Kondo conductance 
in the quantum dot.t2l In the analysis of this experiment 
below, we will see that the frequencies used were of the 
order of Tk, and the range of u> was less than one decade. 
In addition, it was impossible to calibrate the amplitude 
of the field applied to the device; the attenuation coef- 
ficients were frequency-dependent. Therefore, one could 
not perform the measurements of G pca k(w) at a fixed 
field amplitude. However, it was possible to measure the 
dependence of the peak conductance on the amplitude 
of the applied field at a discrete set of fixed frequencies. 
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The ac source was powerful enough to allow the authors 
of Ref. [[2] to completely suppress the Kondo anomaly. 
Nominally, the ac field was applied to the gate, but it 
was apparently hard to exclude, "leaking" of the field to 
other electrodes of the device.E2l It creates further ambi- 
guity in the interpretation of the experimental results. 

The Kondo temperature of the system can be found 
from the dependencies G pea ^(T) and G(Vd c ) measured in 
the absence of irradiation. Comparing the experimental 
results with the theory [Eqs. ( |63|) and fl64|) l , we obtain the 
estimate Tk ~ 30 — 60 /ieV. The external irradiation was 
applied through a high-frequency coaxial cable, coupled 
capacitively to the gate. The frequency of the irradiation 
ranged from 10 GHz to 50 GHz (i.e. Too between 40^eV 
and 200 fieV). 

The zero-bias Kondo peak was clearly observed in the 
G(Vdc) dependence, when no ac field was applied. With 
the increase of the amplitude of the ac field, the height 
of the zero-bias peak decreased for each of the frequen- 
cies used. Such a behavior is in agreement with the our 
conclusion that irradiation must suppress the Kondo ef- 
fect in a quantum dot even when dot ionization does not 
occur. Also, the satellites did not app ear, in agreement 
with our theory for Ticu ~ Tk, see Sec. VIII . 

The authors of Ref. [lj attempted the data collapse 
for the dependence of G pca k on the amplitude V; rr and 
frequency ui of the ac field. In this procedure, Vi rr was 
scaled by some frequency-dependent parameter, in order 
to bring to a single curve the dependencies G pca k vs. V m 
measured at different values of u>. Successful data col- 
lapse means that a universal dependence exists, 



the Kondo conductance by ac bias. This effec t is not so 
sensitive to a specific value of Ed, see Section VII. 



Gpoak = F[Mrr/(M] 
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In the experimental the scaled curves coincided with 
good precision, see Fig. 7 of Ref. [l2| for conductances 
Gpcak ranging from the unitary limit down to values small 
compared with Gjj. Our theory indeed allows for such a 
single-parameter scaling j87|), see Eqs. (70) and ( |7l"| ) and 
Fig. ||, in the case of ac modulation of the gate voltage. 
On the contrary, according to Fig. || scaling of the type 
( [S7| ) would be possible only in the domain of large or 
small conductances, if the ac field of frequency to ~ Tk 
is applied to the leads. 

We must mention here that the experimentally 



measured^ Kondo conductance was suppressed by the 
irradiation uniformly across the Coulomb blockade val- 
ley, including its middle point, where Ed = U — Ed- In 
our theory, however, the suppression is not uniform, and 
vanishes at Ed = U/2, if the ac field is applied to the 
gate [see Eqs. (Q, Q), and (§||)]. We see two possible 
reasons for the discrepancy. First, our conclusion is valid 
for a weak modulation only. This condilipn most prob- 
ably was not satisfied in the experiment^ performed at 
Twj ~ Tk] at such frequencies strong modulation of the 
gate potential, STk ~ Tk, is required to achieve signif- 
icant suppression of the Kondo effect. Second, leakage 
of the ac field to the leads would result in suppression of 



X. CONCLUSION 

We considered the Kondo conductance of a quantum 
dot subjected to ac field. We have shown that, despite the 
essentially non-equilibrium character of the problem, the 
Kondo temperature Tk [Eq. ([37])] remains the only rel- 
evant energy scale. The dc Kondo conductance depends 
on the ac field only through two dimensionless variables 
which are the frequency and the amplitude of the ac per- 
turbation divided by Tk- In terms of these two variables, 
conductance is a universal function. The form of this 
function, and the relation of the perturbation amplitude 
to the "bare" value of the ac field amplitude, depends on 
the way the ac field is applied. 

If the ac field is applied to the gate, then the strength of 
the perturbation is characterized by the amplitude STk 
of adiabatic variations of the Kondo temperature, see 
Eqs. (|5|)-@). At low frequencies, fiu < T K , the con- 
ductance is close to the unitary limit [Sec. |v|, Eq. ([57])]. 
At higher frequencies, hu > Tk, the ac field suppresses 
the Kondo effect by means of the decoherence of the dot's 
spin (Sec. [if], 0). The value of the zero-bias conduc- 
tance decreases with the increase of the frequency u> of 
the ac field. The results we obtained for the modulation 



of the gate voltage are summarized in Sec. VI 



If the ac field is applied across the dot, then the proper 
variable is the corresponding dimensionless bias between 
the leads, V^/Tk- A "slow" field, YuJ < max{eK ac , Tk}, 
suppresses the Kondo effect essentially the same way as 
a finite dc bias does, see Eqs. (76), ([77]). A "fast" ac 



field, hu>' > maxjeVac, Tk }, affects the Kondo conduc- 
tance through the decoherence of the dot spin [Eqs. (|73), 
([73|), (|7^)]. At the fixed amplitude of the field, the sup- 
pression of the Kondo effect diminishes with the increase 
of the ac field frequency to' . 

The ac field also produces satellite peaks in the depen- 
dence of the differential dc conductance on the dc bias. 
However, the satellite maxima in the conductance are 



inevitably small, see Sec. VIII 



The analysis of the experiment^ [Sec. IX demon- 
strates good agreement between our theoretical results 
and the results of the experiment. 
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APPENDIX A: AC BIAS WITH eV» 



In this Appendix we describe the crossover between 
the regimes of "slow" and "fast" ac bias, which occurs at 



(Al) 



see Fig. ||. Throughout this Appendix, we will use the 
finite-order perturbation theory to evaluate the decoher- 
ence rate and conductance. The RG technique is aban- 
doned here, since the finite-order perturbation theory is 
sufficient in the region defined by condition ( |Al[ ). 



1. Decoherence by ac bias 

Unlike the case 7' = eVac/Tiu' -C 1, in the crossover 
region (Al) the decoherence rate is determined also by 
many-photon processes. Using the Fermi Golden Rule 
with the Hamiltonian of Eqs. (|ll]) and (pf), we arrive at 
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where M{x) = -xJ a (x)J 1 (x) + x 2 [J (x)} 2 + x 2 [Ji{x)f . 
In the case of "fast" ac field, 7' = eV^/Tiuj' <C 1, equation 
( [A2| ) reduces to 
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The latter formula is similar to Eq. (41) and accounts for 
single-photon processes only. Equation (A3) is the first 
term in the perturbation theory series in J aa i for Ti/t' . 
The summation of leading terms of all orders in J aa ' 
can be performed with the help of the RG technique, 
and yields Eq. (|7|) [cf. Eqs. (|lf) and 0) respectively] 
In the limiting case of "slow" ac bias, 7' 3> 1, Eq. (|aJ) 
is reduced to 
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cf. Eq. (BO) 



2. Conductance 



The conductance in the crossover region (Al) can be 
evaluated with the third-order perturbation-theory series 
in Jaa>{t). Using Eqs. @, ©, ©, and ©, we ar- 
rive at 
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[cf. Eq. (|24|)1. The expression (A5) accounts for the cur- 
rent induced by both dc and ac biases. To single out the 
former contribution, which is the true dc conductance, 
we must average over the phase cf>o- This averaging in 
the limit of zero temperature and dc bias yields 



G 



(3) 
peak 



2 nh l 



T-(O) 

J LR 



J RR 



7(0) 



iff • 



exp(-|i|/r) 



--Jo 



V* 2 + (H/ttDo) 2 
In the limit w' — > 0, equation ( |A6| ) yields 

GpL « ln(^ /eUac) . 



27' sin 



Lu't 



(A6) 



(A7) 



This result is analogous to Eq. U77J) at eV ac ^ Tr-. At 



small frequencies, oj't < 1, the corrections to Eq. (A7) 



are proportional to exp(— 1/lu't). At larger frequencies 
Ti/t < hio' <C eVac, the two leadin g te rms in the expan- 
sion of the right-hand side of Eq. ( A5 ) in powers of I/7' 
are 



G 



(3) 
peak 



37T 2 e*_ 
32 nh l 



7(0) 

J LR. 



2 r 



J, 



(0) 



RR 



In 



Do 



— lno/r 

7 



At even larger frequencies, Tiu/ ^S> eV ac , Eq. (| 



(A8) 

yields 
(A9) 



This is the first logarithmic term of the series that are 
summed up in Eq. (|7^). 

Note that the results given by Eqs. (|A7j), ( |A8| ), and 
( |A9| ) match each other at the corresponding applicabil- 
ity limits. 



APPENDIX B: RG TRANSFORMATION FOR 
SATELLITE PEAKS 

In this Appendix, we describe the RG transformation 
we used to derive Eq. (|8^). Unlike the other instances 
of application of the RG technique in our paper, the RG 
transformation of this Appendix consists of two stages. 
The first stage is analogous to the one considered in 
Sec. [Ill A| [Eqs. (H|), @]. It stops when the bandwidth 
D reaches eVd c ~ fru. Since S,/r stl < hco [see Eq. (|3l|)], we 



need to reduce the band further to account for all possible 
virtual transitions contributing to the Kondo anomaly in 

Gstl(Vdc)- 

The RG transformation we consider in this Appendix is 
aimed at evaluation of the conductance G st i- In processes 
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which contribute to the singularity in G s ti, a tunneling 
electron jumps from the Fermi level in one lead to the 
Fermi level in the other, emitting or absorbing a pho- 
ton. Therefore, the reduced band in each lead must be 
centered at its Fermi level. When D is below Tiu), transi- 
tions of only two types are possible within such a band. 
First, there can be transitions within a lead without ab- 
sorption/emission of a photon. Second, transitions from 
the higher-potential lead to the lower-potential lead with 
emission of a photon, and reverse transition with absorp- 
tion of a photon are also possible. The other types of 
transitions bring electrons out of the reduced band and 
should be excluded from the consideration. Such a treat- 
ment yields the RG equations 



d^f eta CI y a 

dD ~ V ~D~ ' 

dJhR _ Jlr{Jll + Jrr) 
dD ~ V D 

with the initial conditions 



J aa {D)\ 



2T a 



1 



Jlr(D)\ 



T L +T R 2v\n(tw/T K ) 
1 



R 



ST K 



T L + T R 4 V [ln{huj/T K )}2 X) 



K 



(Bl) 
(B2) 

(B3) 
(B4) 



cf. Eq. (|43|). The second stage of the trans- 
fo rmation must be sto pped at fl ~ D* = 
\J (fi./r s ti) 2 + (eVd c ± fiuj) 2 . Expanding the solution for 
Jlr in powers of STk/Tk up to the first power, we ob- 
tain 



Jlr\ 



d~d* 



2^T L T^5T K 

t l + r R T K 



In 



yW^ti) 2 + (e^dc ± huj) 2 | 



l-K 



K 



(B5) 



The conductance G s t\ must be calculated in the second- 
orde r pe rturbation theory in J^lr(D ~ D* ), g iven by 
Eq. (B5). This calculation finally yields Eq. (§2). 
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